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Peripheral  Milling  of  Thin  Titanium  Plates: 
Modeling,  Analysis,  and  Process  Planning 


Abstract 

This  paper  presents  research  conducted  to  model  and  analyze  the  peripheral  milling  of  thin 
titanium  plates.  Impact  tests  are  performed  and  the  vibration  characteristics  of  a  thin  titanium 
plate  are  modeled  empirically.  The  force  process  is  described  by  a  mechanistic  model  and 
experimental  data  are  used  to  determine  the  model  parameters.  The  Particle  Swarm  Optimization 
technique  is  used  to  determine  the  parameters  of  the  plate  vibration  and  force  process  models, 
which  are  combined  to  create  a  complete  model  of  the  thin  titanium  plate  peripheral  milling 
operation.  The  models  are  validated  experimentally  and  the  results  showed  excellent  agreement 
between  the  simulation  and  experimental  results.  A  process  planning  scheme  for  peripheral 
milling  of  thin  titanium  plates  is  developed  and  applied.  In  this  scheme  integer  numbers  of 
widths-of-cut  and  depths-of-cut  are  used  and  the  optimal  helix  angle  and  feed  are  computed  for 
each  combination  of  width-of-cut  and  depth-of-cut.  The  process  planning  test  case  revealed  the 
optimal  material  removal  rate  decreases  as  the  width-of-cut  decreases,  the  optimal  helix  angle  is 
independent  of  the  width-of-cut,  and  the  optimal  feed  increases  as  the  width-of-cut  decreases. 
The  test  case  also  revealed  the  optimal  material  removal  rate  is  independent  of  depth-of-cut,  the 
optimal  helix  angle  increases  as  the  depth-of-cut  decreases,  and  the  optimal  feed  increases  as  the 
depth-of-cut  decreases. 
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1.  Introduction 


Due  to  its  light  weight  and  high  strength,  titanium  is  a  critical  material  for  a  wide  variety  of 
aerospace  components.  Many  of  these  components  do  not  bear  structural  loads  and  are  being 
designed  to  increasingly  thinner  specifications  to  further  reduce  weight  and  improve  aircraft 
performance.  A  better  understanding  of  the  machining  process  is  required.  Therefore,  this 
research  seeks  to  1)  model  the  plate  vibrations  and  machining  forces  for  peripheral  milling  of 
thin  titanium  components,  2)  validate  the  models  experimentally,  and  3)  develop  a  process 
planning  scheme  for  peripheral  milling  of  thin  titanium  components. 

Literature  Review 

There  have  been  many  studies,  most  of  which  are  experimental  in  nature,  that  have  investigated 
the  machining  of  titanium  alloys.  Malone  [1997]  discussed  machining  tough  materials  including 
titanium  and  suggested  using  inserts  with  high  positive  rake  angles  and  sharp  cutting  edges  or 
light,  precise  hones.  Many  studies,  such  as  that  by  Baeker  et  al.  [2000],  investigated  chip 
formation.  Barry  et  al.  [2001]  used  acoustic  emission  (AE)  to  investigate  the  turning  of  Ti6A14V 
and  found  that  for  high  speeds  the  chip  welded  to  the  tool  and  fracture  was  the  dominant  source 
of  the  AE  signal.  Chen  et  al.  [2002]  investigated  chip  formation  in  high  speed  milling  of  titanium 
alloys  and  found  a  lack  of  chip  defonnation  at  high  cutting  speeds.  There  have  been  several 
studies  of  tool  wear  in  titanium  machining  operations.  Brinksmeier  et  al.  [1999]  investigated  the 
use  of  mist  lubrication  when  end  milling  Ti6A14V  and  found  it  to  be  better  than  dry  machining 
and  as  effective  as  flood  coolant  when  using  high  speed  steel  tools,  de  Lacalle  et  al.  [2000] 
investigated  tool  wear  as  a  function  of  cutting  speed,  feed,  and  tool  coatings  when  machining 
titanium.  They  concluded  that  TiCN  coatings  are  preferable.  Zoya  and  Krishnamurthy  [2000] 
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conducted  end  milling  experiments  of  a  titanium  alloy  using  PCBN  tools  and  observed  that 
700°C  is  a  critical  temperature  that  can  be  a  criterion  for  tool  life,  and  cutting  speed  in  the  range 
of  185  to  220  m/min  is  recommended.  Tool  life  can  be  improved  by  using  liquid  nitrogen  for 
cryogenic  machining  of  titanium  [Hong  et  al.,  200 1  a, b].  Tool  life  was  also  shown  to  improve 
when  milling  titanium  compressor  blades  using  a  cooling  refrigerant  and  a  PCD  tool  [Kuljanic  et 
al.,  1998].  Corduan  et  al.  [2003]  investigated  the  wear  mechanisms  for  machining  titanium 
alloys  when  using  polycrystalline  diamond  (PCD)  tools,  cubic  crystalline  boron  nitride  (CBN) 
tools,  and  TiB2  coated  carbide  tools.  It  was  found  that  the  different  wear  patterns  are  not  only 
strongly  linked  to  the  chemical  reactions,  but  also  to  the  different  phase  changes,  which  occur 
when  using  PCD  and  CBN  tools.  Xu  et  al.  [2004]  used  radial  force  frequency  analysis  to  explore 
a  monitoring  method  for  flank  wear  and  found  that  the  amplitude  of  the  fifth  harmonic  was  very 
sensitive  to  flank  wear.  Wang  et  al.  [2005]  investigated  the  use  of  binderless  CBN  tool  when 
high  speed  machining  Ti6A14V  and  found  that  workpiece  adhesion  and  attrition  are  the  main 
wear  mechanisms.  Another  study  by  Shivpuri  et  al.  [2002]  investigated  the  microstructure, 
particularly  cracking,  of  titanium  parts  after  end  milling  with  a  variety  of  cutting  conditions.  Li  et 
al.  [2004]  used  embedded  thermocouples  to  investigate  cutting  temperature  when  end  milling 
titanium.  They  found  the  cutting  temperature  increased  as  the  cutting  speed  increased.  Li  et  al. 
[2007]  investigated  drilling  operations  for  Ti6A14V  and  found  that  high  throughput  drilling  of 
Ti6A14V  is  possible  when  using  coolant,  spiral  point  drill  geometry,  fine-grained  WC-Co  tool 
material,  and  the  proper  balance  of  feed  and  speed. 

Several  studies  have  modeled  milling  operations  of  thin-walled  structures.  Elbestawi  and 
Sagherian  [1991]  conducted  simulations  to  determine  surface  errors  in  the  milling  of  thin-walled 
sections.  The  end  mill  was  modeled  as  a  cantilevered  beam  and  the  plate  was  simulated  using  the 
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Finite  Element  (FE)  method.  The  study  considered  tool  disengagement.  Tsai  and  Liao  [1999] 
investigated  static  surface  errors  of  thin-walled  parts  in  peripheral  milling  operations  using  the 
FE  method  to  model  both  the  end  mill  and  the  part.  Static  defonnation  of  cantilevered  plates  was 
investigated  by  Aijun  and  Zhanqiang  [2008]  by  considering  the  average  cutting  force  acting  on 
the  plate.  An  analytical  static  deflection  model  was  used  to  simulate  the  plate.  Meshreki  et  al. 
[2008]  created  a  simplified  model,  based  on  Rayleigh’s  energy  method,  to  simulate  milling 
operations  of  thin-walled  structures.  The  FE  method,  coupled  with  impact  tests,  was  used  by 
Seguy  et  al.  [2008]  to  create  a  model  of  milling  operations  for  thin  components.  The  model  was 
used  to  investigate  chatter  and  surface  roughness.  Process  planning  to  regulate  surface  errors 
when  machining  thin-walled  parts  has  also  been  considered.  Budak  and  Altintas  [1995]  changed 
the  tool  immersion  in  the  axial  direction  to  regulate  static  form  errors  of  plates  in  peripheral 
milling  operations.  The  end  mill  was  modeled  as  a  cantilevered  beam  and  the  plate  was  simulated 
using  the  FE  method.  The  effect  of  the  work  piece  tilt  angle  on  the  surface  error  of  thin 
cantilevered  plates  was  investigated  for  traditional  milling  operations  in  Lee  et  al.  [2003]  and  ball 
end  milling  operations  in  Lee  et  al.  [2004].  Zhao  et  al.  [2004]  investigated  controlling  surface 
errors  by  tool  path  optimization  including  tilting  the  tool  and  using  a  combination  of  large  and 
small  axial  and  radial  depths-of-cut.  Ratchev  et  al.  [2005]  investigated  the  use  of  a  tool  profiled 
motion  along  the  depth  direction  to  achieve  a  constant  width  of  removed  material. 

Experimental  System 

The  machine  tool  used  for  the  machining  experiments  conducted  in  this  paper  is  a  Cincinnati  V- 
CNC  750  vertical  machining  center.  The  maximum  spindle  power  is  17  hp  and  the  maximum 
spindle  speed  is  8000  rpm.  The  travel  in  the  x,  y,  and  z  directions  are  762,  508,  and  508  mm, 
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respectively,  and  the  maximum  cutting  feedrate  is  250  mm/s.  A  solid-carbide,  eight-flute  end 
mill  was  utilized  in  all  of  the  machining  experiments.  The  tool  has  a  helix  angle  of  40°,  a  radius 
of  9.525  mm,  and  a  corner  radius  of  3.048  mm. 

To  measure  machining  forces  during  the  operations,  a  3-axis  dynamometer  (Kistler 
9257BA)  was  utilized.  The  ranges  are  ±5  kN  in  the  x  direction  (parallel  to  the  plate  in  the 
horizontal  direction),  ±5  kN  in  the  y  direction  (perpendicular  to  the  plate),  and  ±10  kN  in  the  z 
direction  (parallel  to  the  plate  in  the  vertical  direction).  The  natural  frequencies  for  the  x  and  y 
directions  are  2.0  kHz  and  the  natural  frequency  for  the  z  direction  is  3.5  kHz.  The  amplifier  was 
equipped  with  2.0  kHz  low  pass  filters  for  all  three  directions.  To  measure  plate  vibrations  during 
machining  operations,  eddy  current  displacement  sensors  (Lion  Precision  ECL  100-U5B)  were 
utilized.  The  sensors  are  specifically  calibrated  for  titanium,  have  an  offset  of  0.25  mm,  a  range 
of  1.25  mm,  a  bandwidth  of  80  kHz,  and  a  resolution  of  0.75  pm.  To  determine  cutting  tool  and 
plate  structural  dynamics,  an  impact  hammer  (PCB  Piezotronics  086C03)  was  utilized.  The 
sensitivity  is  2.25  mV/N,  the  range  is  ±2.2  kN,  and  the  -10  dB  frequency  range  is  8  kHz. 

A  National  Instruments  (NI)  PCI-6040E  multifunction  card  was  utilized  to  measure  data 
during  the  machining  experiments  and  for  the  structural  modeling  experiments.  The 
multifunction  card  has  sixteen  single-ended  (eight  double-ended)  analog  inputs.  The  maximum 
input  rate  is  500  kS/s  with  an  operation  range  from  ±0.05  to  ±10  V. 

A  special  fixture  to  machine  thin  plates  of  various  thicknesses  was  constructed  (see 
Figure  1).  The  fixture  is  attached  to  a  dynamometer  for  machining  force  measurements.  Eddy 
current  sensors  are  placed  on  the  side  of  the  fixture,  opposite  to  the  machined  surface,  for  plate 
vibration  monitoring  during  machining  (see  Figure  1). 
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Figure  1:  Fixture  Mounted  to  Kistler  Dynamometer.  Eddy  Current  Sensors  can  be  Placed 
at  Different  Vertical  Positions  with  Increments  of  6.35  mm  and  Five  Different  Horizontal 
Locations.  Machining  Occurs  on  the  Opposite  Side  of  the  Plate. 

Particle  Swarm  Optimization 

Particle  Swarm  Optimization  (PSO),  a  method  based  on  the  concept  of  social  interactions  such  as 
bird  flocking  and  fish  schooling  [Kennedy  and  Eberhart,  1995],  is  used  in  this  work  to  develop 
models  of  the  thin  titanium  peripheral  milling  operation.  Given  measurement  data,  model 
parameters  are  optimized  such  that  the  difference  between  the  model  outputs  and  corresponding 
measurement  data  is  minimized. 

The  algorithm  utilizes  the  following  steps.  First,  the  model  structure  is  determined.  The 
model  contains  np  unknown  model  parameters.  Then,  the  parameters  of  ns  model  parameter  sets 
are  randomly  populated.  Each  model  parameter  is  bounded  by  minimum  and  maximum  values. 
Each  model  set  is  referred  to  as  a  particle.  Next,  changes  in  the  model  parameters  for  each 
parameter  are  randomly  selected.  These  changes  are  also  bounded  by  minimum  and  maximum 
values.  The  cost  function  of  each  of  the  ns  model  parameter  sets  is  evaluated  in  the  next  step.  The 
cost  function  is  a  user  defined  function  that  describes  how  well  the  model  set  fits  the  given 
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experimental  data.  Then,  the  cost  functions  for  the  ns  model  parameter  sets  are  compared  and  the 
model  parameter  set  with  the  minimum  cost  function  is  detennined.  This  model  parameter  set  is 
the  global  best  model  parameter  set.  Also,  the  model  parameters  of  the  global  best  model 
parameter  set  are  the  global  best  model  parameters. 

The  PSO  algorithm  is  implemented  for  a  fixed  number  of  iterations.  The  value  of  the y'th 
model  parameter  in  the  /'lh  model  parameter  set  (xy)  at  the  kth  iteration  is 

xij{k)  =  xij(k~l)  +  vij{k)  i  =  l,-..,ns  7  =  1,...,  np  (1) 

where  Vy  is  the  change  in  the  yth  model  parameter  in  the  zth  model  parameter  set  at  the  kth  iteration 
and  is 


Vj  ( k )  =  w(k) vy  ( k  - 1)  +  eft  (£)[/?..  (k  - 1)  - xy  (k  - 1)] 

+c2r2{k)[pAk-{)-xu{k-1)] 


i  =  \,...,ns  7  =  1,...,  np  (2) 


The  subscripts  i  and  j  refer  to  the  zth  model  parameter  set  and  the  /h  model  parameter, 
respectively.  The  parameter  w(k)  is  known  as  the  inertial  weight,  the  parameters  c\  and  ci  are 
user  selected  constants  (c\  =  C2  =  2  in  this  study),  and  the  parameters  r\(k)  and  ri{k)  are  random 
numbers,  evenly  distributed  between  0  and  1 ,  determined  at  each  iteration.  The  parameter  py  is 
the  best yth  model  parameter  in  the  zth  model  parameter  set  and  p,A,  is  the y'th  model  parameter  in  the 
global  best  model  parameter  set.  Therefore,  the  PSO  algorithm  changes  each  model  parameter 
according  to  the  previous  change  in  the  model  parameter,  the  difference  between  the  model 
parameter  and  the  best  parameter  in  its  model  parameter  set,  and  the  difference  between  the 
model  parameter  and  the  global  best  model  parameter  of  all  of  the  model  parameter  sets.  If  the 
model  parameters  and  changes  in  model  parameters  exceed  their  minimum  or  maximum  values, 
they  are  set  equal  to  their  minimum  or  maximum  values,  respectively 

if  x..(k)<x”m  then  Xj(k)  =  xJm  i  =  l,...,ns  j  =  l,...,np  (3) 
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if  xIJ(k)>x[ r 

then 

r  ffr)-  rmax 
Xij\K)~Xj 

i  =  l,...,ns  j  =  \,...,np 

(4) 

if  vs(*)<vf 

then 

V,  (9  =  vf 

i  =  l,...,ns  j  =  \,...,np 

(5) 

if  v#(q>v“ 

then 

vf(*)=vr 

i  =  l,...,ns  j  =  \,...,np 

(6) 

The  inertial  weight  linearly  decreases  from  its  maximum  (wmax)  to  its  minimum  (rvmm)  value  and 
is 

#)  =  Wma x-kW™*~W™  (7) 

^max 

where  the  parameter  kmax  is  the  maximum  number  of  iterations.  Equations  (l)-(7)  are 
implemented  until  k  =  kmax.  Since  the  PSO  method  is  not  gradient  based,  it  is  capable  of  escaping 
local  minima  and  finding  the  global  minimum. 


2.  Plate  Dynamic  Modeling 

The  structural  response  of  a  thin  titanium  plate  (254  mm  in  length,  50.8  mm  in  height,  and  3.378 
mm  thick)  was  modeled  experimentally.  The  plate  was  mounted  in  the  special  fixture  (see  Figure 
1).  The  plate  was  clamped  on  the  bottom  and  both  sides,  and  was  free  on  the  top.  Ten  impact 
tests  were  conducted  at  fifteen  locations  (see  Figure  2).  For  each  test,  the  eddy  current  sensor 
was  mounted  on  one  side  of  the  plate,  and  the  plate  was  struck  on  the  other  side  with  the  impact 
hammer  directly  opposite  to  the  sensor.  The  impact  hammer  force  and  eddy  current  displacement 
signals  were  sampled  simultaneously  at  a  frequency  of  10  kHz.  Each  data  set  was  windowed 
such  that  500  total  points  (10  before  the  impact  and  490  subsequent  to  the  impact)  were  selected 
for  analysis.  For  each  data  set,  the  displacement  and  force  signals  were  transfonned  into  the 
frequency  domain  using  the  Matlab  fft  command. 
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|  Fixture  □  Titanium  plate  #  Measurement  point:  stiffness  (kN/mm) 


Figure  2:  Titanium  Plate  (3.378  mm  thickness)  Static  Stiffnesses. 


The  magnitude  frequency  responses  at  the  frequencies  of  20,  40,  60,  80,  100,  120,  140, 
160,  180,  and  200  Hz  were  averaged  at  each  location.  These  fifteen  values  were  inverted  to 
determine  the  static  stiffnesses,  in  units  of  kN/mm,  at  each  location  and  the  results  are  shown  in 
Figure  2.  The  magnitude  frequency  response  at  0  Hz  was  not  considered  since  the  values  at  this 
frequency  tend  to  be  distorted  due  to  small  biases  in  the  eddy  current  sensor  signal.  As  expected, 
the  static  stiffness  is  minimal  at  the  top,  middle  location  and  increases  for  locations  closer  to  the 
fixture.  The  static  stiffnesses  are  not  symmetric  about  the  plate’s  vertical  centerline  since  the 
measurement  locations  are  slightly  shifted  to  the  right  side. 

A  model  tailored  to  location  7  is  used  for  subsequent  simulations.  This  location  is  the 
critical  point  since  the  plate  static  stiffness  is  minimal  at  this  location.  The  proposed  model  is 
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where  Di  is  the  plate  deflection  at  location  7,  kp  =  2.78  kN/min,  the  four  identified  natural 
frequencies  for  the  poles  are  960,  1280,  1740,  and  3160  Hz,  and  the  parameters  Q\,  <±2,  £3,  coz  1, 
0x2,  (Ozi,  £p\,  Cp3,  and  Qa,  are  determined  by  applying  PSO.  The  cost  function  is  the  sum  of 
the  square  of  the  differences  between  the  experimental  and  model  magnitudes  for  frequencies 
from  0  to  3580  Hz.  The  data  for  frequencies  above  3580  Hz  are  not  considered  since  there  is 
uncertainty  in  the  frequency  response  near  the  sampling  bandwidth.  The  damping  ratios  are 
bounded  between  10  10  and  0.2,  and  the  natural  frequencies  for  the  zeros  are  bounded  as  follows: 
960  <  a>:  1  <  1280  rad/s,  1280  <  ®z2  <  1740  rad/s,  and  1740  <  a>: 3  <  3160  rad/s.  The  maximum 
changes  in  the  damping  ratios  are  bounded  by  ±0.01  and  the  maximum  changes  in  the  natural 
frequencies  are  bounded  by  ±10  rad/s.  The  algorithm  utilizes  500  model  parameter  sets  and 
executes  1000  iterations.  The  damping  ratios  of  the  best  model  parameter  set  at  each  iteration  are 
shown  in  Figure  3.  These  parameters  remain  constant  after  600  iterations  and  are  £1  =  1.63- 10"", 
£2  =  4.43- 10“2,  £3  =  1.10-103,  Qa  =  8.64-HT3,  fa  =  1.68  10  2,  &  =  6.42-10  3,  and  Qa  = 

_3 

4.21-10  .  The  natural  frequencies  of  the  best  model  parameter  set  at  each  iteration  are  shown  in 
Figure  4.  These  parameters  remain  constant  after  200  iterations  and  are  a>z\  =1190  rad/s,  ®-2  = 
1540  rad/s,  and  coz 3  =  2760  rad/s.  The  natural  frequencies  tend  to  converge  more  quickly  than  the 
damping  ratios  since  the  range  of  the  natural  frequencies  are  bounded  more  tightly.  The  cost 
function  at  each  iteration  is  also  shown  in  Figure  4.  Since  the  model  parameters  are  well- 
bounded,  the  initial  cost  function  was  close  to  the  final  cost  function.  The  cost  function  decreases 
quickly  as  the  natural  frequencies  converge  and  then  decrease  slowly  until  the  damping  ratios 
converge. 
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Figure  3:  Evolution  of  Best  Structural  Model  Damping  Ratios. 
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Figure  4:  Evolution  of  Structural  Model  Best  Cost  Function  (top)  and  Best  Structural 

Model  Natural  Frequencies  (bottom). 

The  experimental  and  model  magnitude  frequency  responses  at  location  7  are  shown  in 
Figure  5.  The  model  magnitude  frequency  response  matches  the  experimental  magnitude 
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frequency  response  very  well  over  all  frequencies,  except  those  near  5000  Hz,  which  is  due  to 
uncertainty  in  the  frequency  response  near  the  sampling  bandwidth.  Figure  6  shows  simulation 
and  experimental  results  for  an  impact  test.  The  results  show  the  model  captures  the  force 
frequency  and  magnitude  very  well.  It  should  be  noted  that  the  fixturing  (i.e.,  specific  clamping 
force  distribution)  will  affect  the  plate  dynamic  model;  therefore,  each  plate  will  have  slightly 
different  dynamic  model. 


Figure  5:  Model  versus  Experimental  Magnitude  Frequency  Response  at  Location  7. 
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time  (ms) 


Figure  6:  Experimental  and  Modeled  Displacements  at  Location  7  for  Impact  Test. 


3.  Force  Process  Modeling 

Machining  experiments  were  conducted  on  a  titanium  block  to  determine  cutting  force  model 
parameters.  Nine  experiments,  listed  in  Table  1,  were  selected  for  machining  force  process 
modeling.  For  all  nine  experiments  the  spindle  speed  and  depth-of-cut  were  Ns  =  1350  rpm  and 
d  =  12.7  mm,  respectively,  the  feed/flute  varied  such  that  0.051  <  f/<  0.102  mm,  and  the  width- 
of-cut  varied  such  that  1.27  <  w  <  2.54  mm.  This  range  is  typical  of  traditional  milling,  as 
opposed  to  high  speed  milling,  of  titanium.  The  experimental  results  for  Case  1  (Ns  =  1350  rpm, 
d  =  12.7  mm ,ff=  0.051  mm,  and  w  =  1.27  mm)  in  Table  1  are  shown  in  Figure  7.  The  figure 
shows  the  forces  acting  on  the  part  in  the  x  (parallel  to  the  plate  in  the  horizontal  direction),  y 
(perpendicular  to  the  plate),  and  z  (parallel  to  the  plate  in  the  vertical  direction)  directions.  The 
forces  are  shown  for  five  spindle  revolutions  in  the  x  direction  and  twenty-three  revolutions  in 
the  y  and  z  directions.  Figure  7  also  shows  the  PSD  for  these  force  signals.  The  time  plots  do  not 
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start  at  time  zero  since  cutting  data  was  presented  only  from  the  middle  of  the  cut  to  show  steady 
cutting  conditions.  In  Figure  7,  there  is  significant  energy  at  22.5  Hz  and  180  Hz,  which  are  the 
tool  rotation  and  flute  passing  frequencies,  respectively,  for  this  experiment.  The  y  and  z  force 
signals  illustrate  the  repeatability  of  the  force  process  during  each  tool  rotation.  The  x  force 
signal  illustrates  the  force  variation  during  each  flute  rotation.  The  variation  in  the  x  force  signal 
from  flute  to  flute  indicates  a  slight  tool  eccentricity  was  present. 


Table  1:  Process  Parameters  and  Mean  Machining  Forces  for  Machining  Experiments 
Used  to  Develop  Machining  Force  Model  (Ns  =  1350  rpm  and  d  =  12.7  mm). 


Case  //  (mm)  w  (mm)  Fx  (N)  Fy  (N)  F,  (N) 


1 

0.051 

1.27 

232 

-291 

-83.0 

2 

0.051 

1.905 

304 

-409 

-118 

3 

0.051 

2.54 

366 

-538 

-156 

4 

0.076 

1.27 

337 

-376 

-127 

5 

0.076 

1.905 

442 

-524 

-181 

6 

0.076 

2.54 

534 

-688 

-246 

7 

0.102 

1.27 

436 

-457 

-173 

8 

0.102 

1.905 

572 

-634 

-250 

9 

0.102 

2.54 

697 

-851 

-344 
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Figure  7:  Machining  Forces  acting  on  Tool  and  Power  Spectral  Densities  of  these  Signals 
(Case  1  in  Table  1:  Ns  =  1350  rpm ,f/=  0.051  mm/llute,  w  =  1.27  mm,  d  =  12.7  mm).  The 
signal  Fx  is  Shown  for  Five  Tool  Revolutions.  The  signals  Fy  and  Fz  are  Shown  for  Twenty- 
Three  Tool  Revolutions. 

To  model  the  end  milling  process  the  end  mill  is  sliced  into  disks  with  height  A z  (see 
Figures  8  and  9).  In  these  studies  up  milling,  as  opposed  to  the  more  conventional  down  milling, 
is  used.  In  up  milling,  the  chip  load  on  the  flute  gradually  increases  from  zero  to  the  full  chip 
load  during  the  flute’s  revolution.  The  helix  angle  is  assumed  to  be  zero  for  each  slice.  Each  disk 
is  analyzed  like  a  face  mill.  The  disk  at  the  bottom  of  the  end  mill  will  be  analyzed  first.  The 
instantaneous  feed  of  the  bottom  disk  of  the  7th  flute  is 

fi  =  f/  cos  (3)  (9) 

where  fj  is  the  feed/flute  and  0i  is  the  angle  of  the  7th  flute  measured  in  the  x-y  plane  from  the 
positive  x  direction  to  the  flute,  with  clockwise  as  the  positive  direction.  The  corresponding  chip 
area  of  the  bottom  disk  of  the  7th  flute  is 

a,  =/,Azcos(6t)  (10) 

The  cutting,  thrust,  and  lateral  forces,  respectively,  acting  on  the  bottom  disk  of  the  7th  flute  are 
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The  forces  acting  on  the  bottom  disk  of  the  zth  flute  in  the  Cartesian  coordinate  system  are 


Fx,  =  -Ft.  eos  ( )  -  Fc  sin  ( 6i ) 

(14) 

Fy,  =  -PTi  sin (&,)  +  FCj  cos ( 6> . ) 

(15) 

tC4- 

i 

ii 

(16) 

The  total  forces  acting  on  the  bottom  disk  are 

Nf 

Fx  =X[-  FTj  cos(6», )+  FCj  sm(0i  )]cr(^ ) 

1=1 

(17) 

Nf 

Fy  =  Z Ft:  sin (0i)~FCl  cos(0l  )}t(0,  ) 

1=1 

(18) 

f, 

(19) 

i=l 


where 


cr 


(4)  = 


o  if  et  <  eex 
o  if  Gen  <  e, 

1  if  e,>e>6„ 


(20) 


The  function  cr(Oi )  detennines  if  the  zlh  flute  of  bottom  disk  is  in  contact  with  the  part.  The  entry 

and  exit  angles  are  detennined  by  the  end  mill  and  part  geometries.  For  the  end  milling  tests 
conducted  in  this  paper  (see  Figure  8),  the  entry  and  exit  angles,  respectively,  are 

*.=§  (2D 
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Figure  8:  End  Milling  Schematics.  Dotted  Line  is  Tool’s  Path  during  Previous  Rotation. 
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The  PSO  method  is  utilized  to  determine  the  force  process  model  coefficients  (i.e.,  Kc, 
Kt,  and  KL).  The  mean  values  for  the  x,  y,  and  z  forces  were  computed  for  the  nine  experiments 
in  Table  1.  For  each  iteration  of  the  PSO  algorithm,  ten  combinations  of  model  parameters  were 
simulated  for  all  nine  cutting  conditions  and  their  values  were  updated  according  to  the  PSO 
algorithm.  The  model  parameters  were  bounded  by  ±10  kN/mirT  and  the  changes  in  model 
parameters  at  each  iteration  were  bounded  by  ±1  kN/mm  .  The  performance  index  is  the  sum  of 
the  squares  of  the  errors  between  the  measured  and  modeled  x,  y,  and  z  mean  forces  for  all  nine 
tests.  Figure  10  shows  the  best  performance  index  for  all  one  hundred  iterations.  Note  the 
performance  index  is  plotted  on  a  base  10  logarithmic  scale.  The  index  decreases  in  a  step  wise 
manner  every  time  the  PSO  algorithm  finds  an  improved  solution.  Figure  1 1  shows  the  best 
parameter  values  at  each  iteration.  The  parameter  values  change  significantly  at  first  and  then 
converge  close  to  their  final  values  after  approximately  sixty  iterations.  The  final  values  are  Kc  = 
1.14  kN/mm2,  Kj  =  2.49  kN/mm2,  and  Kl  =  0.759  kN/mm2.  Typically,  the  cutting  coefficient  is 
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the  largest  of  the  three  model  coefficients;  however,  the  data  is  somewhat  weighted  since  the 
widths-of-cut  are  small.  Figure  12  shows  a  plot  of  the  measured  versus  modeled  jc,  y,  and  z  mean 
forces  for  the  nine  sets  of  data  used  to  create  the  model.  Note  the  solid  line  denotes  perfect 
prediction.  The  average  absolute  values  of  the  percent  errors  for  these  nine  tests  are  12.1%, 
3.27%,  and  3.06%  for  the  x,  y,  and  z  forces,  respectively.  Thirteen  tests  (shown  in  Table  2),  none 
of  which  were  used  to  create  the  model,  are  used  for  validation.  Figure  13  shows  a  plot  of  the 
measured  versus  modeled  x,  y,  and  z  mean  forces  for  the  thirteen  sets  of  data  used  to  validate  the 
model.  Note  the  solid  line  denotes  perfect  prediction.  The  average  absolute  values  of  the  percent 
errors  for  the  thirteen  tests  are  13.6%,  3.00%,  and  5.35%  for  the  x,  y,  and  z  forces,  respectively. 


Figure  10:  Evolution  of  Force  Model  Best  Performance  Index,  Note  Performance  Index  is 
Plotted  on  a  Base  10  Logarithm  Scale. 
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Figure  11:  Evolution  of  Best  Force  Model  Best  Model  Parameters. 


Figure  12:  Measured  versus  Modeled  Mean  Force 
Denotes  Perfect  Prediction. 


for  Data  in  Table  1.  Straight  Line 
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Figure  13:  Measured  versus  Modeled  Mean  Force  for  Data  in  Table  2.  Straight  Line 
Denotes  Perfect  Prediction. 


Table  2:  Process  Parameters  and  Mean  Machining  Forces  for  Machining  Experiments 
Used  to  Validate  Machining  Force  Model. 


Case 

Ns  (rpm) 

ff( mm) 

w  (mm) 

d  (mm) 

Fx  (N) 

F,w 

Fz  (N) 

1 

1550 

0.051 

1.27 

12.7 

226 

-288 

-85.3 

2 

1750 

0.051 

1.27 

12.7 

225 

-293 

-81.8 

3 

1550 

0.051 

1.905 

12.7 

301 

-409 

-120 

4 

1750 

0.051 

1.905 

12.7 

300 

-416 

-120 

5 

1750 

0.076 

1.27 

12.7 

328 

-372 

-123 

6 

1550 

0.076 

1.905 

12.7 

438 

-523 

-193 

7 

1750 

0.076 

1.905 

12.7 

437 

-531 

-189 

8 

1750 

0.076 

2.54 

12.7 

525 

-704 

-260 

9 

1750 

0.102 

1.27 

12.7 

422 

-446 

-171 

10 

1750 

0.102 

1.905 

12.7 

569 

-645 

-261 

11 

1350 

0.076 

1.27 

19.05 

509 

-652 

-220 

12 

1350 

0.102 

1.27 

19.05 

659 

-751 

-295 

13 

1350 

0.102 

1.905 

19.05 

847 

-996 

-420 
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4.  Model  Validation 


The  machining  force  model  developed  in  Section  3  is  used  to  simulate  the  titanium  end  milling 
process  for  Cases  3,  5,  and  13  in  Table  2,  and  results  are  compared  to  the  corresponding 
experimental  data.  The  simulated  and  experimental  machining  forces  for  one  tool  rotation  for 
Cases  3,  5,  and  13  are  shown  in  Figures  14,  15,  and  16,  respectively.  The  model  predicts  the 
shape  of  the  force  profiles  and  the  force  average  values  very  well.  However,  the  experimental 
data  has  some  variation  that  is  due  to  structural  deflections,  which  are  not  included  in  the  model. 
This  is  very  apparent  for  Case  13.  This  case  represented  extreme  cutting  conditions  in  that  the 
feed  and  depth-of-cut  were  at  their  respective  maximum  values  and  the  width-of-cut  was  at  its 
middle  value. 


_ i _ i _ i _ i _ i _ i _ _ _ i_i 
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Figure  14:  Simulated  and  Experimental  Machining  Forces  for  Case  3  in  Table  2:  Ns  =  1550 
rpm,/=  0.051  mm/flute,  w  =  1.905  mm,  and  d  =  12.7  mm. 
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Figure  15:  Simulated  and  Experimental  Machining  Forces  for  Case  5  in  Table  2:  Ns  =  1750 
rpm./  =  0.076  mm/flute,  w  =  1.27  mm,  and  d  =  12.7  mm. 
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Figure  16:  Simulated  and  Experimental  Machining  Forces  for  Case  13  in  Table  2:  Ns  = 
1350  rpm,/=  0.102  mm/flute,  w  =  1.905  mm,  and  d  =  19.05  mm. 


simulation  - experimental 
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The  plate  dynamic  structural  model  for  the  top  middle  location  and  the  machining  force 
process  model  are  used  to  simulate  the  titanium  end  milling  process,  and  the  simulation  and 
experimental  results  are  compared.  The  simulations  are  known  as  Time  Domain  Simulations 
[Tsai  et  ah,  1990;  Smith  and  Tlusty,  1993;  Week  et  ah,  1994].  Results  are  shown  in  Figures  17 
and  18.  These  figures  show,  for  one  spindle  rotation,  the  three  machining  forces  and  the  plate 
structural  displacement  at  the  top,  middle  location  as  the  tool  was  passing  this  location.  In  both 
figures  it  can  be  seen  that  the  simulated  forces  and  displacements  predict  the  general  shape  of  the 
experimental  data  very  well.  For  the  results  shown  in  Figure  17,  the  simulated  x  forces  slightly 
over  predict  the  experimental  data  and,  for  the  results  shown  in  Figure  18,  the  simulated  x  forces 
slightly  under  predict  the  experimental  data.  The  magnitude  of  the  simulated  y  forces  is  very 
similar  to  the  experimental  data  for  both  machining  operations.  The  simulated  z  forces  for  both 
machining  operations  predict  the  average  experimental  data  very  well;  however,  the 
experimental  data  had  substantially  more  variation  than  the  simulated  data.  The  simulated 
displacements  predict  the  variation  in  the  experimental  displacements  very  well  for  the  results 
shown  in  Figure  17  and  slightly  over  predict  the  variations  in  the  experimental  displacements  for 
the  results  shown  in  Figure  18.  The  simulated  displacements  are  also  slightly  larger  than  the 
experimental  data  for  both  machining  operations.  This  deviation  is  likely  due  to  the  fact  that  the 
plate  stiffness  characteristics  slightly  change  every  time  a  new  plate  is  fixtured. 
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Figure  17:  Time  Domain  Simulation  versus  Experimental  Results  for  Machining  Titanium 
3.378  mm  Thick  Plate:  Ns  =  1350  rpm,//  =  0.051  mm/flute,  w  =  1.27  mm,  and  d  =  15.75  mm. 
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Figure  18:  Time  Domain  Simulation  versus  Experimental  Results  for  Machining  Titanium 
3.378  mm  Thick  Plate:  Ns  =  1750  rpm,  /}  =  0.102  mm/flute,  w  =  1.27  mm,  and  d  =  15.75  mm. 
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5.  Peak-to-Peak  y  Force 


The  effect  of  width-of-cut,  depth-of-cut,  and  helix  angle  on  the  peak-to-peak  y  force  is 
investigated  in  this  section.  For  helix  angles  of  5°,  20°,  35°,  and  50°  the  machining  process  is 
simulated,  without  structural  vibrations,  and  the  maximum  of  the  absolute  value  of  the  y  force 
and  the  peak-to-peak  y  force  are  computed  for  a  range  of  widths-of-cut  from  0.254  to  3.048 
mm  and  a  range  of  depths-of-cut  from  2.54  to  30.48  mm.  The  spindle  speed  is  1750  rpm,  the 
tool  radius  is  9.525  mm,  the  number  of  flutes  is  8,  and  the  feed  is  0.051  mm/flute.  For  each  helix 
angle,  three-dimensional  plots  showing  the  maximum  of  the  absolute  value  of  the  y  force  and  the 
peak-to-peak  y  force  as  functions  of  the  width-of-cut  and  the  depth-of-cut  are  constructed.  The 
results  are  shown  in  Figures  19,  20,  21,  and  22  for  helix  angles  of  5°,  20°,  35°,  and  50°, 
respectively.  For  all  helix  angles,  the  maximum  of  the  absolute  value  of  the  y  force 
monotonically  increases  with  increasing  width-of-cut  and  increasing  depth-of-cut.  For  a  helix 
angle  of  5°  the  peak-to-peak y  force  increases  as  the  depth-of-cut  increases;  however,  the  peak- 
to-peak  y  force  first  increases  and  then  decreases  as  the  width-of-cut  increases.  This  trend  with 
respect  to  the  width-of-cut  is  the  same  for  all  four  helix  angles  considered  in  this  study.  For  a 
helix  angle  of  20°  the  peak-to-peak  y  force  increases,  decreases,  and  then  increases  again  as  the 
depth-of-cut  increases.  The  valley  seen  in  the  peak-to-peak  y  force  plot  in  Figure  20  indicates 
the  depth-of-cut  can  be  adjusted  such  that  the  peak-to-peak  y  force  is  very  small.  For  helix 
angles  of  35°  and  50°,  the  peak-to-peak  y  force  plots  have  two  and  four  valleys,  respectively, 
indicating  the  depth-of-cut  can  be  adjusted  to  multiple  values  to  achieve  a  very  small  peak-to- 
peak  v  force.  If  the  peak-to-peak  y  force  plots  were  extended  in  the  width-of-cut  direction,  it 
would  be  seen  that  there  are  valleys  in  this  direction,  too. 
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w<mm)  d  (mm)  w  (mm)  d  (mm) 

Figure  19:  Maximum  Absolute  Value  of  y  Force  (left)  and  Peak-to-Peak  y  Force  (right)  as 
Functions  of  Width-of-Cut  and  Depth-of-Cut:  Ns  =  1750  rpm,  N/=  8,  Rt  =  9.525  mm.  //  = 
0.051  mm/flute,  and  y.  =  5°. 


w  (mm)  d  (mm)  w  (mm)  d  (mm) 

Figure  20:  Maximum  Absolute  Value  of  y  Force  (left)  and  Peak-to-Peak  y  Force  (right)  as 
Functions  of  Width-of-Cut  and  Depth-of-Cut:  Ns  =  1750  rpm,  N/=  8,  Rt  =  9.525  mm.  //  = 
0.051  mm/flute,  and  yr  =  20°. 
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w  (mm)  d  (mm)  w  (mm)  d  (mm) 

Figure  21:  Maximum  Absolute  Value  of  y  Force  (left)  and  Peak-to-Peak  y  Force  (right)  as 
Functions  of  Width-of-Cut  and  Depth-of-Cut:  Ns  =  1750  rpm,  N/=  8,  Rt  =  9.525  mm.  //  = 
0.051  mm/flute,  and  y.  =  35°. 


w  (mm) 
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Figure  22:  Maximum  Absolute  Value  of  y  Force  (left)  and  Peak-to-Peak y  Force  (right)  as 
Functions  of  Width-of-Cut  and  Depth-of-Cut:  Ns  =  1750  rpm,  N/=  8,  Rt  =  9.525  mm.  //  = 
0.051  mm/flute,  and  yr  =  50°. 
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6.  Process  Planning 

A  process  planning  scheme  for  peripheral  milling  of  thin  titanium  parts  is  given  in  this  section. 
The  material  to  remove  has  a  cross  section  with  a  total  depth-of-cut  ( dt )  and  total  width-of-cut 
(w,).  In  this  scheme  the  spindle  speed,  tool  radius,  and  number  of  flutes  is  fixed.  However,  the 
scheme  can  be  expanded  to  include  these  variables  as  well.  The  machining  operation  has  a 
constraint  on  geometric  tolerances  (±0.127  mm)  and  a  desired  material  removal  rate  of  273.1 
mm  /hr.  However,  the  geometric  tolerance  constraint  will  be  considered  by  a  maximum  v  force 
constraint  which  is  given  by 

Fmax  =  2.78  •0.127  =  0.353  lkN  (27) 

where  2.78  is  the  stiffness,  in  kN/mm,  of  the  weakest  point  on  the  plate  and  0.127  is  the 
geometric  tolerance,  in  mm.  This  constraint  will  be  utilized  in  the  process  planning  scheme. 

6.1  Process  Planning  Scheme 

A  maximum  helix  angle  (ymax)  and  a  minimum  helix  angle  (ymin)  are  selected.  In  this  study  ymax  = 
50°  and  ymm  =  5°.  Since  the  number  of  tool  passes  must  be  an  integer  number,  the  process 
planning  scheme  will  consider  combinations  of  depths-of-cut  and  widths-of-cut  that  are, 
respectively,  dt,  dj 2,  <7,/3,  etc.  and  wh  wj. 2,  w/. 3,  etc.  For  a  depth-of-cut/width-of-cut 
combination  the  minimum  helix  angle  is  selected  and  a  bisection  routine  is  used  to  detennine  the 
maximum  feed  such  that  the  maximum  v  force  constraint  is  not  violated.  The  function  used  in  the 
bisection  routine  is  Fmax  -  Fy  m«x  where  Fv,max  is  the  maximum  y  force  determined  by  simulation. 
The  bisection  routine  uses  an  initial  low  feed  of  0  mm/flute  since  this  will  definitely  cause  the 
function  to  be  positive  and  an  initial  high  feed  of  25.4  mm/flute  since  this  will  definitely  cause 
the  function  to  be  negative.  The  algorithm  converges  when  the  following  condition  is  realized 
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F  -F 

3^,  max  m 


F 


<0.01 


(28) 


The  lower  feed  is  selected  to  ensure  the  maximum  y  force  constraint  is  not  violated.  The  helix 
angle  is  incremented  by  A yr  and  the  bisection  routine  is  again  used  to  determine  the  maximum 
feed  such  that  the  maximum  y  force  constraint  is  not  violated.  In  this  study,  A yr  =  3°.  This 
procedure  is  continued  until  the  maximum  helix  angle  is  utilized.  For  the  depth-of-cut/width- 
of-cut  combination,  the  optimal  helix  angle  is  determined  to  be  the  helix  angle  where  the  feed  is 
maximum,  since  this  will  maximize  the  material  removal  rate,  which  is  given  by 

MRR  =  wfdN fNs  (29) 

The  procedure  is  applied  to  every  combination  of  depth-of-cut  and  width-of-cut.  The  solutions 
are  then  checked,  via  time  domain  simulations  for  stability.  The  solution  with  the  maximum 
material  removal  rate  that  is  stable  is  then  selected.  Note  tool  wear  and  surface  finish  constraints 
are  not  considered  in  this  process  planning  scheme.  These  will  be  incorporated  in  future  studies. 
Also,  the  time  to  prepare  for  subsequent  tool  passes  is  not  considered  in  the  material  removal  rate 
computation. 


6.2  Process  Planning  Scheme  Test  Case 

The  process  planning  scheme  is  applied  to  a  test  case  where  dt  =  30.48  mm,  wt  =  3.048  mm,  R,  = 
9.525  mm,  and  Ns  =  1750  rpm.  Up  to  three  passes  in  both  the  depth  and  width  directions  are 
considered.  The  results  are  shown  in  Table  3.  This  example  has  discemable  trends  with  respect 
to  the  optimal  material  removal  rate,  helix  angle,  and  feed.  First,  optimal  material  removal  rate 
decreases  as  the  width-of-cut  decreases,  and  is  insensitive  to  the  depth-of-cut.  Second,  the 
optimal  helix  angle  increases  as  the  depth-of-cut  decreases,  and  is  insensitive  to  the  width-of- 
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cut.  The  optimal  feed  increases  as  the  width-of-cut  and  depth-of-cut  decrease.  In  all  nine  cases, 
chatter  did  not  occur.  The  results  indicate  that,  from  the  perspective  of  material  removal  rate,  the 
total  width-of-cut  should  be  taken.  However,  other  aspects,  particularly  tool  wear  rate  and 
surface  finish,  will  also  need  to  be  taken  into  account.  In  general,  the  lower  the  feed,  the  better 
the  surface  finish,  which  would  indicate  the  total  depth-of-cut  should  be  taken.  Also,  since 
multiple  tool  passes  are  not  necessary,  the  time  between  tool  passes,  which  was  not  taken  into 
account,  would  not  be  required. 


Table  3:  Process  Planning  Scheme  Test  Case  Results. 


Case 

d  (mm) 

w  (mm) 

f  (mm/flute) 

7r  (deg) 

MRR  (mm3/hr) 

Stable 

1 

30.48 

3.048 

1.859- 10“2 

14 

37.59 

Yes 

2 

15.24 

3.048 

3.708- 10“2 

26 

37.59 

Yes 

3 

10.16 

3.048 

5.518T0-2 

38 

37.34 

Yes 

4 

30.48 

1.524 

2.946- 10-2 

14 

29.72 

Yes 

5 

15.24 

1.524 

5.969- 10-2 

26 

29.97 

Yes 

6 

10.16 

1.524 

8.534- 10-2 

38 

28.70 

Yes 

7 

30.48 

1.016 

4.115*1 0-2 

26 

27.69 

Yes 

8 

15.24 

1.016 

8.052- 10-2 

26 

27.18 

Yes 

9 

10.16 

1.016 

0.1179 

35 

26.42 

Yes 

Time  domain  simulations  conducted  for  Cases  1,  3,  4,  and  7,  and  the  y  force,  part 
displacement,  and  PSDs  for  these  two  signals  are  shown  for  these  cases,  respectively,  in  Figures 
23,  24,  25,  and  26.  The  power  spectral  densities  of  the  y  force  and  part  deflection  signals  were 
computed  for  the  last  half  of  the  data  sets,  which  is  during  the  steady  state.  The  power  spectral 
densities  of  the  y  force  signals  contain  significant  energy  at  233  Hz,  which  is  the  flute  passing 
frequency,  and  its  first  and  second  multiples.  This  is  also  true  of  the  power  spectral  densities  of 
the  part  deflection  signals;  however,  there  is  also  significant  energy  at  four  multiples  of  the  flute 
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passing  frequency.  Since  the  power  spectral  densities  only  contained  significant  energy  at  the 
flute  passing  frequency  and  its  multiples,  chatter  was  not  present  during  the  operations.  It  can 
also  be  seen  that  the  feed  is  adjusted  such  that  the  maximum  absolute  value  of  the  y  force  during 
the  steady  state  is  approximately  Fmax  and  the  part  displacement  is  approximately  0.127  mm. 
When  comparing  the  y  forces  in  the  simulations  in  this  section  to  the  y  forces  shown  in  Figures 
14-18,  it  can  be  seen  that  the  peak-to-peak  variation  in  the  y  forces  are  smaller  by  an  order  of 
magnitude.  Therefore,  the  process  planning  scheme  naturally  found  the  valleys  in  the  peak-to- 
peak  v  forces  discussed  in  the  previous  section. 
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Figure  23:  Time  Domain  Simulations  (left)  and  Power  Spectral  Densities  (right)  ofj  Force 
and  Part  Displacement  for  Case  1  in  Table  3:  Ns  =  1750  rpm ,f/=  1.859-10-2  mm/flute,  w  = 
3.048  mm,  d  =  30.48  mm,  yr  =  14°,  and  MRR  =  37.59  mm3/hr. 
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Figure  24:  Time  Domain  Simulations  (left)  and  Power  Spectral  Densities  (right)  of  y  Force 
and  Part  Displacement  for  Case  3  in  Table  3:  Ns  =  1750  rpm ,f/=  5.518-10-2  mm/flute,  w  = 
3.048  mm,  d  =  10.16  mm,  yr  =  38°,  and  MRR  =  37.34  mm3/hr. 


time  (s)  frequency  (Hz) 


Figure  25:  Time  Domain  Simulations  (left)  and  Power  Spectral  Densities  (right)  of.y  Force 
and  Part  Displacement  for  Case  4  in  Table  3:  Ns  =  1750  rpm ,f/=  2.946-10-2  mm/flute,  w  = 
1.524  mm,  d  =  30.48  mm,  yr  =  14°,  and  MRR  =  29.72  mm3/hr. 
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Figure  26:  Time  Domain  Simulations  (left)  and  Power  Spectral  Densities  (right)  of  y  Force 
and  Part  Displacement  for  Case  7  in  Table  3:  Ns  =  1750  rpm,^)  =  4.115-KT2  mm/flute,  w  = 
1.016  mm,  d  =  30.48  mm,  yr  =  26°,  and  MRR  =  27.69  mm3/hr. 


7.  Summary  and  Conclusions 

This  paper  described  work  that  investigated  the  modeling  and  process  planning  of  perferial 
milling  operations  for  thin  titanium  parts.  The  equipment  that  was  used  in  the  experimental 
studies  was  described  in  detail.  The  Particle  Swarm  Optimization  (PSO)  method,  which  was  used 
to  model  structural  displacements  and  machining  forces,  was  described  and  implemented.  The 
structural  models  were  compared  to  experimental  data  in  the  frequency  and  time  domains.  The 
machining  force  models  were  simulated  and  compared  to  experimental  data  generated  by  cutting 
a  titanium  block.  The  machining  force  models  were  combined  with  the  structural  models  to 
simulate  the  machining  of  a  thin  titanium  plate  and  simulation  and  experimental  results  were 
compared.  A  process  planning  scheme  for  the  perferial  milling  of  thin  titanium  parts  was 
developed  and  a  test  case  was  conducted. 
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The  PSO  method  was  successful  in  creating  machining  force  models  and  structural 
vibration  models.  The  force  models  predicted  the  machining  forces  very  well  when  machining 
thin  plates.  The  structural  models  had  some  variation  from  the  experimental  results.  This  is 
believed  to  be  due  to  the  fact  that  the  plate  lixturing,  which  is  slightly  different  each  time  a  plate 
is  put  into  the  fixture,  affects  the  structural  models.  The  simulations  showed  that  variations  in  the 
y  force  during  a  flute  rotation  were  insensitive  to  width-of-cut,  but  are  very  sensitive  to  depth- 
of-cut  and  that  for  specific  depths-of-cut  the  variation  was  virtually  zero.  The  process  planning 
test  case  revealed  the  optimal  material  removal  rate  decreases  as  the  width-of-cut  decreases,  the 
optimal  helix  angle  is  independent  of  the  width-of-cut,  and  the  optimal  feed  increases  as  the 
width-of-cut  decreases.  The  test  case  also  revealed  the  optimal  material  removal  rate  is 
independent  of  depth-of-cut,  the  optimal  helix  angle  increases  as  the  depth-of-cut  decreases, 
and  the  optimal  feed  increases  depth-of-cut  decreases. 
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